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ABSTRACT 

The field of Peo electronics has become one of burgeoning 
importance in recent years due to the vast strides made in the field of 
communications and other weTctod fields. Theories abound, but few 
impregnable mathematical descriptions of the phenomenon of the positive 
column have been advanced. An attempt was made to extend one of” the 
mathematical solutions to a more goneral case. The resulting equations 
were programmed into a digital computer which produced results in- 
dicating the efficacy of the program, but more data will have to be 
digested in order to establish the reliability of the basic equations. 

The writers wish to express their erate appreciation to 
Professor Norman L. Oleson for suggesting this problem and for encour- 
acing us and guiding our research. We wish to thank Professor H. 
Martinez for his guidance and instruction on the computer portion of 
this thesis--without his help this thesis could not have been brought to 
a successful conclusion. We would also like to thank our two most 
invaluable assistants--our wives, whose unflagging cheerfulness, faith, 


and understanding made this work possible. 
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CHAPTER I 
INTRO DUCT ION 

This thesis evolved from the analysis of a paper by Professors Me 
Satosi Watanabe and Norman Le Oleson of the U. S. Naval Postgraduate 
Schoo1 ‘?? in which the phenomenon of traveling density waves in a 
positive column was Peoria ted: It seems appropriate in the beginning 
to define exactly what is meant by some of the torms and expressions that 
will be used heroaftere This is done expressly to avoid confusion over 
semantics and not as a gratuitous collective insult to subsequent readers. 

ite ke simplest instance let us assume that we have a cylindrical 
glass tubo filled with a pure rare gas at a pressure of say, 12 mm Hg. 
In terms of every day experience this would bo considered a relatively 
low pressure; however, most authorities in tho field of gaseous conductors 
refer to this as a high pressure. Sealed in either end of our gas tube 
let us imagine that we have two electrodes which conform to the shape of 
the tube. Let each of these electrodes be connected to an external source 
of eemef. through an external resistor. Now, it is well established that 
if either the eem-f. or the resistor is varied, a glow discharge can be 
induced to occur in the tubee Generally speaking, the glow discharge is 
a steady state gascous discharge which appears atreduced pressures and is 
charactcrized by relatively low current densities and relatively high 
potential gradients; perhaps it is better to say in this latter case, 
relatively high values of - electric field to pressure ratios, about 
Which more will be said later. This gaseous discharge has a typical 


structure which is marked by a rather steep potential gradient at the 








cathode and operates basically by electron liberation resulting from 
positive ion bombardment at the cathode. Bike simply, this amounts to 
nothing more than having the gas act as a conductor. However, the pheno-= 
mena occurring within the gas while it is acting as a conductor cannot be 
so simply explained. In fact so many phenomena occur that it is not 
foasible to discuss all of them ae the scope of this thesis. As 
indicated in the opening statement, we shall confine our attention to a 
single phenomenon in the positive column of the glow discharge. In order 
to facilitate our development we shall make some simplifying assumptions 
about the composition of our positive column. We shall assume that it is 
composed of throe fluids: 
a. A negative ion "gas" which we shall further assume to be 
composed entirely of electrons | be 
be A positive ion "gas" which we shall further assume to be 
composed either of atoms or molecules of the gas with 
which the tube is filled, but each of which carries a 
single positive charge 
ce The neutral particles of the gas itself--these particles 
are oy er in the majority. 
It is customary in discussing gaseous conductors to refer to ionized 
regions containing high and approximately equal concentrations of positive 


(2) 


and negative ions as the plasma With the proper environment there can 
be caused to appear in the positive column alternate dark and light bands. 


These bands may move or they may be stationary. in either case they are 


called striations. If they appear to move, they are called moving 
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striations; if they do not move, they are called standing striations. 
As a parenthetical mattor of intcrest, it seems well to note that moving 
striations are not detectablo by the unaided eyee 

The traveling density waves referred to earlier may or may not be 
moving striations. Certainly ono could be excused if in the first glow 
of understanding, he impetuously ascribed the title of moving striations 
to these traveling density waves; however, sober and mature consideration 
would make such a conclusion doubtful. The important point to remember 
here is that moving striations are an observable phenomena and not an 
ethereal speculation. 

A mathematical approach by Watanabe and Oleson has been undertaken. 
i their approach they attempt to provide a theoretical basis for ex- 
plaining moving striations. However, they assert that certain of their 
approximations may be in error if applied to actual Peiations >. It is 
the intent hereafter to justify or refute these assumptions. In order to 


do this it is proposed to extend the linear treatment of Watanabe and 


Oleson to the non-linear case and thereby arrive at a more general solution. 
\ 


A search of the literature was instituted to garner the latest, most exact 
data available. Where these data were found to be sither doubtful or non-~ 
existent, data wore taken from current experiments being performed at the 
Ue S.- Naval Postgraduate School. The ultimate aim was to substitute these 
data in the pertinent mathematical equations and solve them by means of a 
digital computer. The solutions then could be compared with the results 
obtained experimentally and a reasonable conclusion could be reached as to 
the validity of the original premises. 
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Since we have assumed that we have a fluid filled tube, let us estab= 
lish a coordinate system to describe the particles in our fluid. in this 
system the neutral particles are always fixed in space and only the positive 
and negative particles move with respeot to our coordinate system. 

if @ is the density of any physical quantity, Q the rate per unit 
volume at which the quantity is being produced, and §V the volume occupied 


by the quantity; then, 


QSV 
or Se paw giv 


Let us define $V as a small volume of fluid and agree that $V always 
designates a volume element that moves with the fluid sush that it contains 
the same fluid particles. in general the volume, 8V, will then change with 
time. If we assume that $V is in the form of a reotangular box, then let 
us say that its dimensions are Sx, Jy, and $z, such that; 

§V= Sx by Sz 
Now, the x=component of fluid velocity, vx, may be different for the op-= 
posing faces of our box moving in the x-direction. If so, 4x will change 


with time at a rate equal to the difference between these two velocities; 


thusly, cae 7 VV F 
joey °7 


and by a similar development we can say, 


a Mi d ~ oy 
Sy = and az $2 = sp Sz 


~{: 





itereiore, 


ae 
aS Sy $2 BS) + 6x $2 % Sy + Snby & 2 


= (2% 4 8%, av 
(Set att U2) sx by Se 


Gee ve el 


f 


Inserting this in the equation above for ait we have: 
= eC a 
SV se + ev VIV=QV 


dividing out 4V, 
dP loy.7- 
me? may v= @ 


but by definition, 


se ot = 
1 rea G7 


therefore, 
c ~ a 
a rv: VEO VEG 
Syne 
Which by a relation in vector analysis reduces to: 


ac = 


This of course is nothing more than the equation of continuity. A more 


: : : : : ** 
rigoreus derivation can usually be found in most texts on mechanicsSe 


, : 
See for instance K. R. Symon, Mechanics (Addison-Wesley Publishing 
Company, Inc., Cambridge, Massachusetts, 1953), p. 265 ff. 








For the sake of simplicity we shall assume that all of the negative 
ions are electrons and that the positive ions are formed by the collision 
of electrons with neutral particles and by no other process--a reasonable 
assumption if we know that our gas is not an electro-negative eae see 


fige 1 for a simple graphical representation of this ionization phenomenone 


> 0 | Oe 


Fig f (Towi2eo Gasas, vow Eweer) 


If we place a charged particle in an electric field, the particle will 
move under the influence of the Coulomb law of attraction. It seems 
reasonable that if there are neutral particles in the field which are un= 
affected by Coulomb forces, the charged particles moving in the field will 
collide with the neutral ones. This obviously ocours. However, since the 
charged particle is still under the influence of the field, the particle 
will again tend to move after collision toward the electrode opposite in 
sign to the particle. The overall effect is to give the charged particle 


— 


a motion generally in one direction (see fig. 2). 


o+°..6 ere boa 
: oe 


Fig. £ ( (Iows2 60 GASES, VOW Eee) 


' Since the Coulomb forces operate to move the particle, it would 


appear that the same law acting between oppositely charged particles would 
cause them to recombinee This doos happen under proper conditions; ieee, 
high pressures (of the order of one atmosphere). However, at low pressures 
the collision frequency and collision losses are low and the relative 
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velocity of the charged particles is higher. Because the velocities are 
higher, collision without recombination eas place since the attracting 
field is not powerful enough to deflect the ions and make them move 
closely together over a large enough Twn 

Let us now define: | 


n = negative ion density 


nf = positive ion density 


<j 
it 


negative ion drift velocity 

or = positive ion drift velocity 

Z = Piriter of lon pairs created per unit time per one ; 
electron (constant if the temperature variations 

(1) ) 


are very small 


then from the equation of continuity: 


gn nv) - 2020 

ot pe 

eqn ‘Oy yn’ +9: (n' V+) -21- O 
ot 


Further consideration of a charged particle in an electric field 
will inevitably lead to the conclusion that the velocity of the ion is a 
function of the field. Experimental evidence shows that the velocity is 


proportional to the field to a first approximation, viz: 


V~E 
(see figs. 3 and 4) 
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Figure 3. Drift velocity of electrons in argon as a function of E/p. 
Re A. Nielsen, Phys. Rev. 50, 950 (1936) 
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Figure 4. The drift velocity of Argon getomic ions in Argon as a function 
of B/p « ds A. Hornbeck, Phys. Rev. 84, 615 (1951) 
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The use of vector notation is justified when consideration is given to the 
fact that the motion of the ion is directed. We may therefore say: 
Thaw 7KE 

where K is a proportionality constant called the mobility--a quantity which 
will be more fully discussed later. Considering again figs. 3 and 4, it 
can be seen that the magnitude of the velocity of the electron is much 
greater than that of the positive ion; indeed, in certain instances the 
velocity of the electron can be shown to be a thousand times as great as 
that of the positive ione This fascinating facet of the properties will 
be given more light later in the discussion of ambipolar diffusion. 

Envisioning a cloud of ions (or electrons) in the tube, one is led to 
the inevitable conclusion that the charged particles may exist in varying 
densities. If this be the case, then obviously the ions will flow from 
regions of higher density to the regions of lower density by diffusion. 
This will of course contribute a component to the velocify at which the ions 


move, and we can say: 


—_ 


Veteeae ai) ~~ Vn 


where Vn is the concentration gradient with the negative sign indicating 
that the motion occurs in the direction of decreasing concentration. 


Therefore, 


ed 


em — D 
Voireusiow ~~ 1 vn 


where D is a proportionality constant called the diffusion coefficient for 
the gas. A more adequate treatment of the derivation of this last 


equation can be found in some detail in most texts dealing with the subject 


) 





of diffusion. Combining these contributions to velocjty we have, 
~ De 
eae y= KE > Vi 


The electric field under which the particles act is composed of two 
parts, the impressed eleotric field and the electric field generated by 


the space charges; such that, 


F=£,+€’ 


‘ 


where E, is the impressed eleotric field and will be considered constant, 
and Et is the electric field generated by the space charges and depends on 
the number of electrons and ions present; i.ee., the density. Therefore 


we can say (from Poisson's Equation): 


VE=V-(E,tE')= WE’ = 47el(n-n) 


where e is the electronic charge. 

It is proposed to solve equations (1) by assuming that the actual 
solution is a superposition of a steady uniform solution n# plus a non~ 
steady, nonuniform disturbance 3). In order to linearize the equation 
Watanable and Oleson assumed that the nonsteady nonuniform disturbance was 
small compared to the uniform, steady disturbance, which is nae in agree= 
ment with experimental aoe (6) In the following development no 
such assumption will be made--to the contrary, we shall assume that at 


some particular time the orders of magnitude of the steady, uniform and 


*See for instance L. Be Loeb, Basic Processes of Gaseous Slectronics 
(University of California Press, Berkeley and Los Angeles, 1955), p.188 ff. 
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nonsteady, nonuniform parts may be approximately equal. This is not to 
imply that they are always equal, but rather that they could be equal at 


some particular time. The word "uniform" as used here means that the 
quantities do not depend on the x-coordinate taken along the axis of the 


tube (1), In keeping with our proposed solutions let us now write: 


~ 


ee): 
eqn (3) ye 
a A a 2 
vie define E€ as the field generated by the steady, uniform charge 
densities ng. 


Substituting the value given for V in equation (2) into equation (1): 


as V' [nr (-KE = [z vn} -2nN=0 
ont sO [or(Kte - {Font} - 200 


clearing the parentheses: 


20°49. nKE-(D 9 })- 2nz0 
20+ [nt KE ~ {p' yn} ]-2n=0 


Substituting equation (3) in the preceding equations: 


ae : eee eee. eo 

set 22 4u. En, KE-w KE-Don- Dvn J- 207-27 =0 | 
gn,” sa es tote nt Ventas tf y= “-Zp =0 

Set SPTVINKE tVKE Di yni-D'y rh f-2N, 


Directing our attention to ‘the terms containing n> and recalling that by 


definition these are constant; therefore, 


IN, t i: 
ee and D VNno=O0 


a 








which leaves, 


See 7 En Ke - UKE-D Va ]-2n- 2220 


Q/7 


== 


we 6 gCnfk'E UKE Don ]-2n;- 22520 


but we have stipulated that E = Ey + E'; it then follows: 

yp ee er SS  apet ae ee 

Spt V EM KEM KE’ KE-~UKE-D Vw ]-2n,-Zu 20 

22s Ve [ntK'E, tek Et eK E,* 2K ED] -2n.- 2 =0 

Again we refer to the fact that Ey 1s a constant, consoquently, 

VEQ=0 

Omitting the zero terms, our equations become: 

= tO [-n KE -DKE,-v Ke'-D vu [- 2n,- 2:0 

a Ve [nine tet KE, tL K'E'-D' Vn ]- Zn, -2V =0 


We have defined € as the field generated by the uniform, steady charge 


density aye Let us now define € as the field due to the nonuniforn, 


nonsteady disturbance y*; such that, 


E's €+€ 


inserting this relation into the equations above: 


oS En Kem KE- 2 KE.~ UKE-DKE-D VP ]-2n,- 20=0 


——_— 


OU Lat KE tnt KEV KE, 1 DK E ty E- DVL] - 2Me- 20-0 
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Removing brackets, 

I 

TE NK VE - KV: B- KE Oz - KV VE-KE-VU -Ki7.e-Ke.7y7 
~0 V*L-2n>-22°=0 

Ser a CDE ante C+K" ey iy +  f toe 

t v%e KV CHKE VY tKVOOCTKE Vo eKyig CrK eV)" 
See) ty ae 

From Poisson's Lquation: 


Vic =47€(n-V, ) and Vez 4mre(v-y) 


where e is the charge on the electron. Recalling that our region of 
interest is near the axis of the tube removed from the ends and assuming 
cylindrical symmetry, we find: 


t oe 
Vpt = 2s and €=€, 


With these conditions in mind let us rewrite the equations: 

a2 nt K [4are (nbn; Jc Lh we (o zt KE, e ~ Kai [Yre (n-n7)]- KE 
ai yay [4 7e (v?- p)I- K€, -D a 

Ipt + o ~ t _2pee 

so tne K Lu melny-ne)]t ae tenetit-y +KE, 32 2" +K'v"[A we (ns-n,)) +K ESB 


rk pt [are (view) tke, 2a # Ye - zn-227=0 


Inasmuch as our interest in the field is in the central region along the 
axis of the tube, we note that the field € (or more particularly €, ) will 
be zero because of the steady, uniformity of the charge density ae causing 


is approximately equal to n_. 


ite Further, in this region n’ ¥ 


O 


gb = = a hi, =O 
-D2=-KE YD - ¥p— "-D)iyr- 2 


“| 
x | 


eqn (4) 
=] -p* 32 K'E, 2 apy “YU rk G, 2 32" ey (Biz) Zn, =O 


Ld 








wnere , 


y24meKk (2n,-N) )-z ¥,247eK N+ 2 
i “ste Kn, Y; = 4ire K 

Tt 
Y,= 4reX” (2ni-1) 2 47eK 


ror comparison purposes the terms, 

~K @, 32 -Y, (2-27) ar - BN, 

tKE, au Y (V-2)2*- any 
which appear in equations (4) do not appear in the original paper by 
watanabe and Oleson. These, then, are the terms dropped by Watanabe and 
Oleson as a result of the mathematical expediency employed by them to 
linearize their equations, and as a result of disregarding certain higher 
order terms. (1) 

Equations (4) are the ones for which constants will be computed for 


a particular set of conditions; eeg., Argon gas, 12 mm Hg pressure, etc. 
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CHAPTER IT 
EVALUATION OF CONSTANTS 
Reiterating our original contention, we proposo to apply equations (4) 
to an experimental situation. Im order to do this certain constant terms 
must be evaluated. Other constants are well knovm but will be repeated 


here in order that our results may be identically reproduced. In this 


; 


latter category of constants we find: : 
e = 4.8 X 10710 CeSeUe 
oees3 X 105-° erea/ °K 


f= 3 1416 


In the former class of constants there will be a need to lmow the values of 


Kv oer pv Plan 24 Dn, ond nv. Rather than list them arbitrarily, it 
seems more appropriate to derive them individually in order to indicate their 
reliability. 
The particular experimental situation to which we shall apply equations 
(4) is as follows: 
A tube 0.5 cm in radius containing pure Argon at a 
pressure of l2é mm Hg with currents of from 20ma to 


(6). 


100ma in our gaseous conductor 
With these conditions the value of the applied electric field oan be found. 


From £1g6 5 it can be seen that for a pressure=radius product of 6 we have 
0 

in our case an a ratio of approximately 0.5. At a pressure of 12 mm Ilg 

Eisenia yield ai ik, of 6 volts/cm. The value of BQ which we obtained for 


=2 


our case yields a value of 2X10” statvolts when converted to e.seu. 


* . 
For a description of the method of obtaining these curves refer to the 
source cited in connection with figure 5. 


LS 








From probe measuroments in the plasma within the tube it has been 


experimentally determined that (6). 


ae ue s 7x 1010 ions/om® 





(cm mm Ho) 


Eee ¢ 
(Ionized Gases, von Engel) 


It appears rather obvious that a relation exists between mobility, 
diffusion, and the number of ion pairs created per unit time per one 
electrone Since the latter two terms have been defined, it remains only 
to define the first torm--mobility. We have rather arbitrarily stated 
that mobility is the proportionality constant which relates the drift 
velocity to the field. However, it seems more important to have a 
"feeling" for the phenomenon of mobility than to have a pedantic lmowledge 


of the word itself. Consider azain a charged particle moving under the 
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influence of a field and within a region occupied by the neutral particles 
of a (Pigs 2). We oan say, thon, that mobility is the drift velocity of 
an ion moving in a unit field.” 

In order to utilize the experimental data available (§), it is neces- 
sary to determine values of mobility for the positive ions and the 
electronse In dealing with positive ion mobility it was found that there 
is a difference in mobilities depending on whether the ions present are 


"atomic™ ions or "molecular" ions. It has been shown that when pressures 


of the order of 6 mm Hg are encountered in Argon, AS and at ions are found 


*x 
In searching for a best value for mobility the phrase "reduced 
mobility” was often encountered. Reduced mobility is best described 
as\/): 


275 


k= ti(Teol 


bIi< 


where p is the gas pressure and T is the gas temperature in °K. Many 
results are given in terms of reduced pressure Doe (see for example 
fige 4). This reduced pressure is defined as (8)$ 


Very often in studying graphically presented results one observes 
that experimenters and theoreticians will plot ion mobility versus 
field to pressure ratio (i.e., ©) -«where the pressure may be reduced 
pressure or that of the gas; P eese, figs. 3 and 4 It can be 
shown(2) that p is inversely proportional to the mean free path, 


thusly: 


ee 


r\. 
If this be so, then: : 


E 
E EA 


where EA is an expression for the energy acquired in a mefepe 
This gives the researcher a more vivid and acsurate picture of what 
is happeninge 


iy | 








to be present in the ratio of 0.5 (5), The mobilities of these two ions 

wore found to be 26 and 1.6 Sone respootively\?), Assuming that the 
volt/om 

two types of ions are distributed in the ratio indicated throughout the 


positive column, a weighted average mobility was found to be: 


feo 10+ om?/stvesec 
If the two types of ions are assumed to be present in equal concentrations, 
then: 
K7 = 4.38 X 10# om’/stvesec 
As can be seen, the orders of magnitude are the same or at Icast near 
enough that either torm could be used with some desree of accuracy. 
From fige 5 the drift velocity of an eleotron was determined to be 


0.56 X 10° om/sec. To a first approximation, 


5 CeO 


Koce 
ee lie ar 


imate joe slice om°/stveseo 


Cl I< 


The diffusion coefficients, pt and D, are most easily found by 
employing the mobility. The relationship between mobility and diffusion 


is derived in numerous texts on the subject (2), The relation is: 


mh oan 
kT 


. 
D 
where oe is the electronic charge, k is Boltzmann's constant, and T is the 
ion temporature in °K. The temperature of the ions is approximately equal 


(10), From this equation, then, the 


to the temporature of the gas 
following values were determined: 
pv = 5048 om” /sec 
D” = 6.46 X 10* om*/sec 


In order to compute the value of z it is necessary to employ ambie- 
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polar diffusion in the positive colum results from the previously shown 
fact that the diffusion coefficient for electrons is greater than for ions. 
Assume that as an original condition there exists near the axis of the 
positive coluwm an almost equal number of electrons and positive ions. 
After a short interval the electrons will tend to diffuse toward the walls 
of the tube at a rate greater than the positive ions. This will cause a 
negativo charge to be set up on the walls of the tube. The preponderance 
of positive ions romaining near the axis of the tube (duo to their slower 
diffusion) causes an electric field to be established between the axis of 
the tube and the wall. This field in turn will cause a slowing down of the 
olectrons moving toward the wall and a speoding up of the ions. It happens 
that wall recombination is preponderant at the pressure we have specified. 
This is borne out by the absence of recombination spectra indicating that 
volume recombination is not important here. In an equilibrium state the 
losses due to recombination are equal to the ratos of production. The 
Picoeric field across the positive column is such that electrons arriving 
from the faraday dark space aro accelerated to produce enough electrons 
and ions per unit length of the tube to compensate for losses, thus 
establishing eetd ibrium <7. Specifically, in our case we have restricted 
our consideration to the region near the axis of the tube so that redial 
variations neod not be considered. The point to be remembered here is 

that under the action of the radial electric field, an averace velocity 


of diffusion is soon reached such that, 


t ang cng’ 


Vv a v —— = 
Yr Yave dr dr 


ros” 
which leads to the fact that tho total ambipolar current is zeroe As we 


have alreadyseen, 


ng = no 19 





Also, we have previously implied (page 13 ) that near the axis of the 
tube 6, =» O and €, = 0. Using the expressions for velocity and eliminating 


the field between them, it can be shown iaaeya a). 

pv me px? 
. aS 
where D, is defined as the ambipolar diffusion coefficient. Further, 
since KX is much greater than reiere I, 1s much greater than T. (where dips 
and T, are the electron and ion temperatures respectively), it can be 


shown that: 


T 
DMT DY Gyag D1 x7 - Kt 
I» K™ 8 


Employing knovm data (6) (a= desoc0% and Ty = 300 °K) in the approxi- 
mate formula for D, yields: 
= 145 om”/sec 


From the ambipolar solution to our problem it can be shown that: 


ae eda 
(R)* 


where R is the tube radius. Inserting the value for D, computcd above 
end the radius of our tube=-0,.5 ere for the conditions we have chosen, 
we find: 

2 = 3.35 X 10° ion pairs/electron-sec 
Inserting the constants just computed into equations (4) yields: 


ow 


a —£.44 X10 sa = 2X10" <2 2a +761 KIO'(y- py") = 19 X10" e, 22 - 095 (of-v) i ~2,348 x00 


yx 


" 3.4 Bs Bos 3 Se XIO(L" -23) + 4.04 xjo" Ge £2.44 X10" (pt-p) Dt 2 ays xiblt = 


Ae) 





Little accuracy is sacrificed by disregerdinzs those valuos of constants 
which are dominated by the exceedingly largo values of other constants 
which are additive. 

In order to solve this rather formidable differential equation, 
little error is involved if the partial derivatives aro approximated by 
total derivatives and minutely small increments are taken to be used in 
a digital computer. This may be more clearly seen later where the above 
equations are set up for programming into the computer. Once the equations 
are written in a "useful form" for the computer, we find we must have some 
starting point; i.se., some initial values fori, ce ex.» ctce at some 
time, te. Since the approach we have taken appears to be a rather unique 
one, there were no readily available precedents by which we could be 
guidede Harking back to the orizinal intent of this thesis, the answer 
suggested itself. Since Watanabe and Oleson ‘1) assumed a solution to 
their linear equation, it seemed reasonable to take their evaluation and 
apply it to our problem. If this satisfied our equations, then credence 
is lent we the original hypothesis thereby extending the theoretical con- 
tention to actual conditionse Watanabe and aieeon assumed linear 
solutions: 

. oe ave ikx-ivt 

w)ts ave Lkx=iwt 
where a ane a are complex amplitudes, k is the wave number, andw is the 
angular frequency. Thus, another problem arose--the feeding of a complex 
number into a digital computer. It is not impossible to execute such a 
program, but in view of the physical construction of the computer, it would 


be necessary to sacrifice some of the accuracy of the oricsinal provram in 


a 


order to do so. Howevor, since our solution consists of a real part and 
an imaginary part; that is to say, a real solution and an imaszinary 
solution, either is a solution to our probleme So, we will program only 
the real solution. This reduces our solution to: 

yp - a cos(kx-ut ) 

sjt = attos Clox=it ) 
Letting the time arbitrarily be zero initially these reduce to: 

DM # a cos kX, 

vie ee kXo 
Again we have introduced terms which must be evaluated. Directing our 
Penton to k, the wave number, and referring to the work of Watanabe 
and aipeen / we observe that we must assume their linear counterparts 
of our equations (4). This results in no inconsistencies when the 
restrictions placed on our equations are considered. Watanabe and Oleson's 


equations are: 


SU eye) O20) \ ee - . 
28 Ke, (32) - 9 BL + yar- ywtze 


t xe 


l= D' 22! + ¥,V"- ae 


where t, Y;, y,, and Y, are the same as the ones previously defined. 


Substituting 
2 = ae Lkx= Air 
i) = avo us 


into Watanabe and Oleson's equations and performing the indicated 
operations will give: 
(-iveikK-B k f DK? AY, Jam = Y,a7 = 0 
-yan f (iw ik78 k + pi? J, at = 0 
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In order for two homogeneous equations to have a solution the determinant 


of coefficients must be zero. 


-iw-ikK EQk ADK” of Y, - ¥, 
=40 
— ¥y -iwe iK7S ole f Deed Y; 
Upon solving this determinant we have, 
[-iw- ik Bok AD Ke HY Yi] | -iw Aik? Ek J oh ¥3] - ¥2%4 = 0 


Recalling that if a complex equation can be written: 
RA# iIl=0 
where R is the real part and iI the imaginary part, then R = O and iI = 0 
for the equation to be satisfied. 
Hore particularly for our case, 
horn =O, 


Ke AKER kee UKE gk fk KEP a (Dk eof %) = ¥, 8y 


| 
SP, 


for I = 0 (since i # 0), 
sopphPd i] age (oh? £4] ofr? 44] HitoeLor? 44] = o 


Solving this last equatidn forw: 


iy [e2(K/or - KDA - Ty - TY, I 
fot Arte (4 4%) 


Eliminating G@) between the two equations yields, 


eqn (6) (D7Ké / ¥, \(pt 2 a y, )- VU, 
1 BB ee la ts ore? £h) - 
[om 4G) x (oe? 4 4]? 
Introducing new variables, 
eqn (7) e: Dk’ Af, and H = on? Y; 
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Solving for k* in the first of equations (7), 
{ 
22 h-hh 
~ 


Substituting this value for k* in the second equation (7), 


pf 
aa = (§ - t,) 
Taking cognizance of o fact that D™ » i can say: 


(me. 


Referring to the original expression for n ; 


a - ph.2 $ y, 
and considering the preceding approximation, we observe: 
pr? 2 0 
Since pt is not zero (see page 18), this approximation can be true if and 
ovguly gsr 
k® 2 0 
Inserting the approximate values of 7 and ¢ from equation (7), et seq, 


into equation (6) and observing that k* 2 Gp 


pe Shey 


if we map equation (6) in the yn plane, we have approximated very closely 
a point on the extremely complicated curve represented by equation (6); 


(1) 


2 
1eGo, the point where k = 0. It can be shown from an elementary 


geometrical relationship that kmowing a point and the slope of a curve at 


the point, the equation of the tangent can be written. In our case, 
(I+ 
eqn (8) fa |" 6) (S-g} 
ee agra 8 i) 


a4 





= \2..2 
where, Q&= baly and h = xt oa ake 
YS D 


Substituting the values of and : as shown in equations (7) for the 
appropriate terms in equation (8), and noting that: 

aK 47 eK*n 
we solve for ae ‘ 


Kt 


x7 po wr pr of KZESKE 
7 7 4ime 






Recalling the proposed linear solution, we direct our attention to the 


complex amplitudes and write: 


at = roid (with A and g real) 


a 


which can be written: 


ae ~ A(cos @ - isin @) 


& 


but as has been shown, 
(iW = iKTEQk Jo“ + Ve Y, Ae = 0 
From these last two equations we can say, 


7 é - i: a ee 
at - Acos gb » jAsin g = a 
a 2 


Equating reals and imaginaries, 


ee , Wa) 
hoos f= DE z and Asin gs 
; | 





Solving for A, 


ae Th en ae ee rae J tan°o 
Y, cos to 
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rrom which, 
() # K7Eok 
ml es, 


4) = ee! (xt egesian ak 2. 
[(o* AD ee 4 #4,)] 


Simplifying yields, 


tan £ = 


Inserting, 


(xt wr K™ )kEo 
ean (9 tan @ = 
an (9) Pri ere yO FD 


It can be shown ‘1) that tan # has one maximum between O and Vee 
Taking the derivative of tan J with respect to k and equating to zero we 
find at this point: 

[2b0(K7 ~ K) | (pt Y Ie 


Squaring equation (9) and inserting equation (10) gives: 


eqn (10) k 


(x7 San 
[aoe # OY, x YY 


Calculating this value of @ from the know values of the constants 


(tang) = 


in the preceding equation, we used this as the maximum value of a phase 
factor with which to enter the computer along with our previously described 
boundary conditions. The value of # obtained as described was found to 


be 1°54", 
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CHAPTisR III 


ASSUMPTIONS AND CRITIQUE OF ASSUMPTIONS-- 
SOME DIFFICULTIES eNCOUNTERED 


The experimental data yield a Langmuir plot which is a straight line 
over tvo oycles leading to an assumed Maxwellian distribution for electrons \°) 
The assumption that the ions have a Iaxwellian distribution stems from the 
fact that at several millimetcrs pressure, the tomperature of the ions is 


approximately equal to that of the gas molecules 10), 


The assumption that 
we have radial ambipolar diffusion undor our stated conditions within the 
tube is based on tho results of Karge, Hooks, and Oleson‘$), In this work 
they state that the instantaneous field intensity is due to a very small 
surplus of positive ions or electrons. Their reasoning for this is very 
emply stated in their work.- Since we oan say that n & woes it follows 
that we can say the ambipolar current at the walls is approximately zoro. 
Further, 


dn 


ae — 
Poe oo. = Y 


ei 5 
where we define V+ and vt as tho negative (iee., eloctron) ion and 


positive ion radial drift velocities respectively. 
In our preceding work and the subsequent computer work it was assumed 


that zg was a constant. If zis also evaluated from the yolaticd cue 


1 
eqn (11) z2£9X10 ape Emme V; 


(6) 


using the temperatures observed by Karge, Hooks, and Oleson* ‘, it is 


found that z varies widely with temperature as shovm bolows 


LG z 
16, 440 1.89 X 10° 
14,660 5.3 X 10° 
8, 200 eee 


al 





(11) notes a variation in z based on temperatures (the data 


Farris 
being taken from Pupp's work). Vhen this information is applied by the 


uso of our equation (11), the following results are obtained: 


Me z 
28, 000 : Daal ol 
11,520 es 6 hee 


it should be noted that the validity of the equation for z depends on 
having liaxwellian distribution of elootrons which we have assumed. 


Farris (1!) 


has recently published a short paper on a theory of moving 
striations based on the diffusion theory. His basic assumption is quoted: 

If it is now assumed that the small difference between the 

ion and electron concentrations has the same effect axially 

as radially, and that mean values of mobilities and diffusion 

coefficionts may be used, the theory of a stationary positive 

column may be readily gencralized to include changes of 

concentration axially and with time. climinating the current 

end ficld from the usual conservation and current density 

equation gives: 5 s dn/dt = 0. én, whero S is the rate of 

production of ions per unit volume per second and D, the 

Schottky ambipolar diffusion coefficient given by: 
The major assumption different from our assumption is that the axial 
drift velocitics of the clectrons and positive ions are assumed equal 
within the positive colum. Using our equations (1) and oliminating the 
field wo obtain: n¥ = D,Vne Substituting this in our equation (1) we get 
Farris! oquation (1). Ilowever, we question the validity of the assumption 
that the axial drift velocitics of all ions in the positivo column are 
equale The intcrestins point is the variation of z with electron 
temperature o 


rrom these considerations (®) we observe that there is a fallacy in 


our assumption that z is a constant. Truo, we imposed the condition that 
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this was valid only for small temperature variations. But at the time 
we imagined this variation could easily be on tho ordor of 100 °K. Now, 
it appears that this is not so. Certainly this is a field for further 
inquirye How this deviation from our original assumption affects our 
problem will be shown in the next chapter. We have taken zng to be a 
constant of the valuc of 2.345 x 104, 

It seems noteworthy to point out that thero is a phase difference 
between the maximum degree of ionization and the maximum dezree of 
excitation which has beon obsorved in some rare easos, - however this was 
not observed by Karge, Hooks, and Oleson (6), AS B@pein. in fact. thoy 
observed that a maximum in light intensity was accompanied by a maximum 
To, and a minimum in light intensity was accompanied by a minimum Tg. 

This would suggest that thero is no phase difforence in the case of Argon 
ions in Argon gas at 12 mm Hg prossure.e If we can assume that equation (11) 
holds and further assume that the maximum light intonsity is also evidence 
for maximum excitation, then we can say indeed there is no phaso differ- 
ence between the phenomena. 

The temperature dependonce of mobilities and diffusion coefficients 
was ignored in this thesis. Since the evaluated constants in equation (4) 
were taken to be average “hares h it is suspected that temperature 
dependence is of little significance. 

The relationship between mobility and diffusion usod in this thesis, 


e , (pagel8) may be a tenuous one. Von pmeen? states that the 
kT 


oe 
—_ 


ss 
D 


* 

Private communication with Ne L. Oleson, who with Ae We Cooper, has 
found that maximum light intensity end maximum electron temperature are 
not coincident in neon at several mm pressuro. 
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tho oquation is not applicable to electrons for at least two reasons; 
electrons drift with velocities which are not proportional to the applied 
fiold and their mean free path depends on thoir speede Im defense of our 
use of this equation we relied on experimontal data (see fig. 3) which 
shows that for our conditions the velocity is approximately proportional 
to the appliod field. The E obtained by Karge, Ilooks, and Oleson was 
6 volts por om per mm lige P hig figure is well within the range of ficlds 
wherein exporimental data reveals that velocity is proportional to ficld 
(fige 3). No attempt has been made to justify or refute the mefepe 
difficulty. Von Enge1 4) further states that the equation holds in the 
case of ions which move and diffuse in moderately strong fiolds: in 
strong fiolds the velocity is proportional to tho square root of the field. 
Ile furthor states that it can be assumed that the equation holds for weak 
fieldse That our conditions fall within the purview of Von Engel's state-= 
ments can be seen from figure 4, wherein the velocity of the ions is shown 
to be proportional to the field. 

Graham and Ruhlig (12) obtained a value of 2 X 10° om”/sec for the 
diffusion coefficient of the eloctron (D") by using Droyvesteyn distri- 
bution in conjunction with the Boltzmann transport equation. This value 
was calculated at the same conditions specified for our casee As can be 
seen from page 18 we calculatcd a value of 6.46 X 10% cm/sec. No 
attempt was mado to resolve this difference. Ilowever, as will be showmm 
later, this difference may bo rather casily inveszizated by use of the 


f 


program we propose by merely varying the coefficients of certain torms. 
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CHAPTSR IV 


PROGRAMING 


In order to facilitate an understanding of the computer program we 


will set down a brief mathematical description of what is taking place 


within the computer. For convenience let pa z and ee W--equations (4) 


may then be written: 


pha E K7(B, fe,)-23 -{[% 7 ee -w)] z 


/Y,w # 2.345 x 10*4 


ray AK(E, # 0x) 3= - | ¥ - 47 eK" (2 - w) | W 


x 


¢ ¥,2 2.345 X 10°* 


and from Poisson's ign, 


bY 
ee 
gx 


= 47 e(z - w) 


Furthermore, we arbitrarily establish as an initial condition: 


as 


Let us now divide 


will apply into n 


Let : 


Let: 
Zi 


Lula 


Cy. 


- QO, at x a 0 
the part of our positive column wherein our equations 


equal parts. 


denote the origin of our column ( = 0) 


be the first point at a distance Ax from x, 
= Lie Wert e041, «ose, te 


= Xj - X3_1 2 L/n, L being the length of the colum 


l= 


= AGE Ale ) 
= w(x4,t ) 


= 6x(x;) » 47%e(2z; - wy) (from iii.) 


ol 








Remembering that it was proposed to replace the partial derivatives with 


total derivatives 


mye une finite appr 





Utilizing this exp 


we have: 

az. 
: a 

16 G:) an 

dw. 

dt 

Leys e,, 
at 


It is proposed to 


we restriot our pr 


Rearranging equati 


dz. 


(page 21 ), we now observe the technique of doing this 


oximation formulae: 


wy. = U. 
LAL aL 


al Ax 


iA be a 
x; ) 


(ax)* 


it 
I}e 


ediency on our first three equations in this chapter, 


wm iE 
= ear ye 45717 224 Fa el) = FB Fen, (Ai fr - 23) 


-[Y; f 4 eK7 (25 - wy)fz, # ¥y wy 420345 X lig 


= ——s (Ww, fy - on #Ws 4) ¢o (0 fo) - Ox Xe, eye ¥,) 


(Ax )* 


2 oe le 
-[¥, - 40K Ge w, Jw; 4 Y, 2; #24545 X 10 


=o, | A aHre(ax) (25 = m4) 


L@= 
examine our solution over ten increments; therefore, 
oblem as follows: 


pen Oasis fie eee 


W112 211 = 9 


ons isa), ii.a), and iii.a) above, we have: 


ee 
| oe * ERM 7 °x3)| ; 


02 








+ (A, tex.) - ae - YH anek? (2, - wo] Zs 


(Ax) 
aw = 
ii.b) ~ = Fe. EG #ox,)| ¥i x1 


7 Fea 2501 7 Uy wy £26345 X 107° 


(ax)* AX 


-{Ee, Z Ox) Lees + LY, ~ 4ieK (2, - ")] Ww; 
aires ieee f Y,2, 42.345 X 1074 


iis Sa en ae 4q9e(ax)(24 = wy) 


Let us now introduce some new constants by combining some of those which 
frequently occur; specifically let, 


A 0587 (arbitrarily setting x g 3cm) 


1 = 


am 10 
Ef 
| 


ho = KZ 2 1.35 x 104 


AneK? = 2.44 X 1074 


0 


Az 


DD 
By s ———~g = +720 X 10* 
(4x) 


pe 3 ae! 


Bo = 


! 
> 
va 


Bz = 4ffek” = 109 


Inserting these new constants into equations ib), iieb), and iiieb) 
yields: 

Payer 1A, = ADE, Hex.) 

ae patos 2 Oy Fens sizes 


A {halo H0n,) = thy ~ Aslig # (ar = mI] ay 


Aiyz321 7% Agniw, 72.345 x 10*4 


Oo 





iio) = > (3, - BoE, 4 ox wi fa 
-[39(8, ~ ex.) #23, + Bgl mp ~ (2; = ve )] fe, 


# Bywi-y pf Bengzi 2.345 X 10/4 


-9 
iii.c) Ox, = Ox; 1 6.023 X 10 “(ax)(2, = wy) 


It should be noted that the value of ae is arbitrarily takon to be 7 x Were 
While this may or may not be exact, it is the correct order of magnitude. 
Further, in the ambipolar region ae s oe « i order to avoid overcomplicating 
the manipulation of decimals a new variable is introduced: 


Za = 1015; 


; Ws 104%, 


Changing variables and inserting the values of the constants where applicable: 
dz. 5 sm ye Ia : , om / 
= {,387-,135 x10" (02+ Ox) (2. + 9.195 X10 (024 Cx,) ~, 77¢ —, 244 X10-L7 Kio + (zim) z 
dt ) 
+, 387 2,_, +,1708 X10* We +, 2345 


‘ 


Ww 
wees = 


7 / ‘ 
r [.720 x10", 6x10" (.02 + Ex;)]w!, +f x10" o2rex;) + 1.44 X108+.109 X10" 7x10" (2, gy I} Wi 


+ 720X10" Wy +, 763X102" +, 2345 


nae ee 475° / 
Cx;> Cus, £1807 K108(z: - ws ) 


OF 





t t 
% us define “a = 002 we 7 and oF = 2; - Wi, then tho preceding 


equations can be written as: 


dz. 
i — Cad 5 ’ ; (yaaa) —_ a ; 
7 = [389 -,135 X10 5. 125,,¢ £135 K10°5)7~ 1174 = 244 X10" [7 X10 5,5 PZ’ 


4,387.2" + 1708 X108 w. +, 2HS 


Ae ' 
ye =[,72Z0K/0"—, X10" Sy Jw, t £6410"); + LAP KIO r 109 K10" (7 X10-1- $; I} W 


+, 720 K10" wy, +762 X10°2) +, 2345 


Cy, = Oxy + 1907 X 10% 52; 


These are the final forms of the equations which the computer will be 
asked to solvee A goneral idea of the method of solution can be obtained 
by considering these equations in conjunction with the flow diagram which 
appears on the following pagee The details of execution of the program 


are found following the flow diagram. 


ers) 





Cx;70 
j— 
1807 X 10° 83; 
Cx;_, 1807 X 10" $2; = Cy; 
-OfL t€z 
3 Se KIO: conn 
339-735 x 105S,; 
&, 2 (.389-.135 x/055,;) 2%, 
ies Srp 
-. 244 X10'7(. 7X10" + Sai) 
~174 - (244X102 (.7 x10" + S23) 
135 x10°S,; - 774 +(.244 2107 (7 XI0°' t Sas) }f 
Q2=(.135x10°S,¢— C74 Let x10? (. 7x07 +50; J} Bye, 
a, 7.387 x 2;., 
Ay =,1708 ¥ 10° w; 
a,t dz 
A3+(a,+ 32) 
Ayt[a, +la.raz)] 
Seay +[azt(a,+az))}} 
mS eiek Si 
,720 X10%= .b KIO" Sj 
b, = (6.720 x 10%- . 6X 107 5,3) wh; 
17 x1074 = $93 
109 x10 (,7X10-4 = S,: 
144 xio%e (.109 x 10°7(,7 x107% = $2; )] 
6x 107 Sst {144 X10" (109 X10 (.7x10-4 = S2i)] 
by= wi x [:6.x/07S)i + {1.44 K10%% [109 X/0'%(. 7x10" "= Sz) )}] 
b> =.720 X10" xX Wy., 
by= - 763 x/0° 2 
b+ b, 
b,+(b, + 6,) 
by +{b, + (b, + bz) 
or ,2345+{b,+[b, +(b,+6,)] 
{ri —) 
rest mt! >} 
ho 


yes 


Numerical Integration Sub- routine (Runge- Kutta- Gill method) 
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The following is the detailed program of the computer solution to our 


problom. 
Loo Opn Opn M BM hi Remarks 
Code A é S 
0010 md 26 2100 2100 0205 0-6, 
O0Oll ad 35 2100 0202 0200 1x22 41 x e244 fix atey i x 3724 


0012 ad 35 0200 0200 0207 

0013 ad 35 0210 0207 0014 

ool4 su 36 [ 24 Wy 0202| 24 - Wy = 8,4 in 0202 
0015 tm 34 0202 2100 0020 no +0016, yes ~0020 
0016 md 26 2100 2100 2000 0 S001 


e200L 
0017 tm 34 3000 2100 0022 


2000 


OGZo sf 31 0240 0202 2000 3 
e807 Aalto Sgj in ee 


CozE mr 25 O24 2001 2001 


0022 ad 35 2100 0205 2002 
- ence 
e Xjo] 2005 
0023 ad 35 2100 0206 2003 
| .. (2000 
0024 tm 34 3000 2100 0400 gives e,, in esi 


.02 in Ee 


0025 ad 35 2100 0242 2002 
2003 


0026 ad 35 2100 0245 2005 


Bes 


0027 tm 54 S000 2100 0400 woe e.. «= S,. in 
f Ay Li 2001 


B15 in 0203 and 0204 
0031 ad 35 2100 2001 0204 


0030 ad 35 2100 2000 mn 
0032 tm 34 2001 2100 0035 no 70035, yes ~ 0035 
0033 md 26 2100 2100 2000 O-— 2000, 2001 
0034 tm 34 3000 2100 0037 
0035 ad 35 2000 0244 5 

~.135 X 10 $43 


0036 mr 25 2001 0245 2001 


oT 





0037 
0040 
0041 
0042 


0043 


0045 
0046 
0047 
0050 
O0S1 
0052 
0053 
0054 
0055 
0056 
0057 
0060 
0061 
0062 
00635 
0064. 
0065 
0066 
0067 
0070 


OO72 


ad 


ad 


md 


ex 


ad 


ex 


md 


ad 


ad 


ad 


ad 


ad 


ad 


md 


ad 


ad 


ad 


09 
O95 
04 
04 
26 
04 
O2 


05 


54 
26 
04 
o4 
OO 
O65 
05 
O98 
O5 
05 
o4 
o4 
26 
04 
35 
20 
39 


55 


2100 
2100 
3000 
2001 
2100 
3000 
0014 
0212 
2004 
A] 
2100 
39000 
3000 
2100 
2100 
2100 
2100 
2100 
2100 
39000 
2000 
2100 
3000 
0252 
0255 
2100. 


2100 


0246 
0247 
2100 
2100 
2100 
2100 
O2ul 
0215 
O2tt 
2100 
2100 
2100 
2100 
2000 
2001 
2100 
0202 
0250 
O251 
2100 
2100 
2100 
2100 
2000 
2001 
0254 


0255 


2002 
2005 
0400 
0045 
2000 
0054 
0212 
2004 
0050 
0053 
2000 
0054 
0500 
0214 


0215 


2000 
2001 


a, 
2003 


0400 


Q066 


2000 


0072 


2000 


2001 


2002 


2005 


58 


e989 in 2002 and 20035 


2389 - .135 X 10° sy; in 4 


no ~0043, yes+0045 


no-O0051, yes+0055 


-. (2000 
i Be 


. (0214 
ae te 


S84; in 2000 and 2001 


7 X 107* in 2002 and 2003 


07 X 10°* #84; in 2000 and 2001 


no-+0064, yes ~ 0066 


pe enel ena % LOT" Ate, 1) 


e774 in 2002 and 2005 





0072 
0073 
0074 
0075 
0076 
0077 
0100 
0101 
0102 
0105 
0104 
0105 
0106 
0107 
0110 
Onli 
0112 
0113 
0114 
0115 
0116 
0117 
0120 
OLl2l 
0122 


0125 


6X 


md 


ad 


ad 


ad 


ad 


su 


eX 


md * 


ad 


ad 


ad 


Ox 


md 


o4 


OL 


o4 


26 


O4 


oO4 


OO 


59 


od 


39 


o6 


O2 


of 


O4 


O4 


Ts) 


o5 


OL 


o4 


26 


of 


o4 


3000 


Oelg 


[ 21 J 


2100 


3000 


0000 


2100 


2100 


2100 


2100 


2004 


2004 


[41-1] 


2100 


0000 


0000 


2100 


2100 


0212 


2004 


2100 
3000 


0000 


2100 
0211 
2100 
2100 
2100 


2100 


2000 
2001 
0257 
0212 
0215 
0211 
2100 
2100 
2100 


2100 


2100 
2001 
0226 


G21 


[ wy | 2100 


2100 


2100 


2100 


0400 


0074 


0077 


2000 


0101 


0510 


at 
0217 


2001 
2004 
2004 
0107 
OPlZ 
2000 
0114 


0516 


0220 
0221 
2004 
0120 
0125 
2000 
O125 


0524. 


=el74 =.244 % J0r* (.7 X 107* ¥ 3.4) im 
2000 and 2001 


no--0075, yes — 0077 


0216 
ag in bez 


e087 in 2001 


no +0110, yes~0112 


az in 0220 and 0221 


az in 0220 and 0221 


no~0Ol121, yes +0125 


og 





0124 
0125 
0126 
0127 
0150 
01351 
0132 
01335 
0154 
0135 
0156 
0157 
0140 
0141 
0142 
01435 
0144 
0145 
0146 
0147 
0150 
Olol 
0152 
0153 
0154 


0155 


ad 
ad 
ad 


ad 


ad 


ad 


ad 


ad 


ad 


ad 


ad 


ta 


pr 
ht 


OX 


ad 


ad 


39 
39 
35 
39 
34 
OO 
Ss 
o4 
39 


39 


30 
35 
o4 
35 
o4 
33 
ol 
22 
30 
30 


Of 


39 


oF 


2100 
2100 
2100 
2100 
5000 
2100 
2100 
3000 
2100 
2100 
3000 
2100 
2100 
3000 
2000 
3000 
O157 


2000 


2001 
0014 
2004 
0154 
2002 


3000 


2000 
2001 
0220 
0221 
2100 
0216 
0217 
2100 


0214 


2100 
0227 
2100 
2000 


0147 


Fy 


2000 


0227 


0156 


0277 


2100 


2100 


0222 
a, in 0222, 0223 
0223 
2002 
az in 2002, 2003 
2003 
0400 az # a4 in 2000, 2001 
2002 
ag in 2002, 2003 
2003 
0400 ao a3 Aaa in 2000, 2001 
2002 
a, in 2002, 2003 
2003 
0400 ay Aaof%az fay in 2000, 2001 
2002 
02345 im 2002, 2003 
2003 
0400 d2i sn 2000, 2001 
au s 
2000 dzi dz; een 
ar nS scaled) in 2 
0145 
0150 
0002 
e 
dz; 
AS + in 2002 


a 
2004 ae in mz of 2004 
0154 

0154s fz} | in mg of 0122 


transferred to intesration sub-routine 


40 





OL96 
0157 
O160 
O161 
0162 
0165 
0164 
0165 
0166 
Or67 
0170 
O171 
O1L7Z 
0173 
O1L74 
OL75 
O1L76 
O177 
0200 
0201 
0202 
0205 
0204 
0209 
0206 


0207 


md 


ad 


su 


ad 


ad 


md 


ad 


Ox 


o4 
26 
54 
55 
29 
36 
55 
55 
o4 
o4 
26 
o4 


35 


o4 


5000 
0203 
0204 
2100 
2100 
2100 
5000 
2001 
2100 
5000 
0120 
2004 


3000 


OOO] 


2100 
2100 
2100 
0264 
0260 
2001 
0266 
0267 
2100 
2100 
2100 
2100 
Onis 
0211 


2100 


O00] 


7777 
0001 
0163 
2000 
0170 
2000 
2001 
2001 
2002 
2003 
0400 
0174 
2000 
OS11 
2004 
05 04: 


0504 


no 70161, yes +0163 


il 
6 X 10 844 


; in 2000, 2001 


6720 X 10* in 2002, 2003 


.720 X 10* ~.6 X 10” sy; in 2000, 2001 


no-~0172, yes ~0174 


1X a7te £4 x 2744 


Sj tomp storage 


j Sale 
ts 


eee os foo) 


41 





0210 su $6 1505 1306 0202 dummy for 0014 


0211 Tt r if extracted for 0054 
0212 mr 25 f e200t- 2001 dumny for 0036 
0213 0002 1e re 


0214 
&, temporary storage 
0215 


O216 


ao temporary storage 
0217 


as 
0220 
az temporary storage 
u) 


0224 
0222 

a, temporary storage 
0225 
0224 ; 

dzi temporary storase 
0225 ao 
0226 0001 s if 
0227 02 if f, 0050 
0250 

by storage 
02351 
0252 

bo storage 
0253 
0254 

bz storage 
02355 
0256 

by storage 
0237 
0240 i f 00535 

+1807 X 108 
0241 1047 3476 0000 


42 





0242 
0245 
0244. 
0245 
0246 
0247 
0200 
0261 
0252 
0255 
0254 
0255 
0286 
0257 
0260 
0261 
0262 
0265 
0264 
0265 
0266 
0267 
0270 
OCF. 
0272 


0275 


02 


02 


02 


Ol2l 


S227 


v071 


QO000 


0441 


6142 


3061 


ep al 


1700 


2670 


0000 


1604 


0062 


0407 


7500 


OT77 


2601 


2226 


0967 


2939 


1156 


1750 


ke 


4061 


6600 


0000 


Q000 


one 


Hee 


t 
202 
0000 
OOLT7 m 
e135 4 10 
7720 
f 
} — 
0142 
fr 
! 7 X 1074 
3173 
0047 
; =2.44 X 1022 
6333 
f 
=~ 114 
1362 
f 
0387 
4671 
0033 - 
Bulyio) ake 
0000 
f 
0 2345 
1156 
0030 : 
ene 
0000 
0017 
~720 X 10% 
0000 


.109 x 1019 


1644 X 10+ 


4.5 





0274 
0275 
0276 
0277 
03500 


O501 


0521 
0322 
0525 
0324 
0325 


0326 


md 


ad 


ad 


su 


ad 


ad 


ad 


md 


ad 


ad 


ad 


04 


26 


04 


o4 


o6 


36 


58 


38 


39 


04 


o4 


26 


o4 


35 


29 


ers) 


35 


i 


£ 


0041 


7066 2107 7440 


0015 


ig 


2100 
5000 


5000 


2100 
2100 
2100 
2100 
2100 
2100 
5000 
2001 
2500 
5000 
2000 
2001 
2100 


2100 


00135 


sg 


[wi] 2100 


2100 
2100 


2100 


2000 
2001 
0202 
2100 
0250 
0251 
2100 
2100 
2100 
2100 
0270 
0271 
0272 


0275 


3 


0027 


0307 
2000 
O311 


0540 


0250 
0231 
2001 
2000 
2002 
2005 
0400 
05235 
2000 
03550 
2000 
2001 
2002 


2003 


.763 X 102° 


m—~1l ml f, m= 12 (oceal) 


(Qn41) xX 2796 


no — 03505, yes 70507 


by in 0230, 0231 


siete = an 2002, 2003 


o7 X 107* = soi in 2000, 2001 


no--S21, yes ~3523 


slog) x 102" x (.7 X 107* = 844) 


1044 X 10* in 2002, 2003 


4A. 





0527 
035350 
0331 
0552 


05535 


0540 


0541 


0542 


03435 


0544 


0554 
0355 
0556 
0357 


0360 


tm. 


md 


ad 


; 


ad 


ad 


ad 


ad 


ou 


Ox 


ad 


ad 


o4 
o4 


26 


O90 
29 
o4 
Of 
o4 
26 
o4 


o4 


ors) 


o5 


09 


05 


o6 


52 


O4 


26 


o4 


o4 


O90 


39 


5000 
0204 
2100 
9000 


0205 


2100 
E2100 
2100 
2100 
0120 


2004 


2100 
9000 


9000 


2100 


2100 


2100 
2100 
2100 
2100 
0264 
0265 
2100 
O21l 
2100 
2100 
2100 


2100 


2000 
2001 
0266 
0267 
0215 


O211 


[¥i-1 | 2100 


2100 
2100 


2100 


2000 


2001 


0400 
03535 
2002 
0335 
2002 
20035 
0400 
0337 
0542 
2000 
0544 


0550 


0252 
0255 
2000 
2001 
2004 
0552 
0355 
2000 
0557 


0560 


0254 


0255 


1.44 X 10° .109 X 101" xX (.7 X 10°*-sgyy 
in 2000, 20001 
no~0331, yes +0333 


no~+0540, yes-~0542 


bp in 0232, 0233 


»720 X 10* in 2000, 2001 


no +0553, yes 03555 


bz in 0234, 0235 


45 





O3561 


0563 
03 64 
0365 
0366 
0567 
0370 
0571 
0372 
03735 
0374 
0379 
0576 
0577 
0400 
0401 
0402 
04035 
04.04. 
0405 
0406 
0407 
0410 
O411 
0412 


ad 


ad 


6x 


md 


ad 
ad 
ad 


ad 


sl 


ad 


ad 
to 


sf 


su 


59 
o& 
Of 
o4 
26 
o4 


54 


55 
55 
35 
o4 


o4 


rae | 
58 
4 
o6 
50 
55 
O7 


ol 


56 
30 


2100 


2100 


0074 


0274 


0275 


0211 


[ 24 | 2100 


2100 
5000 


5000 


2100 
2100 
2100 
2100 
3000 


5000 


3000 
0425 
5000 
2002 
2008 
20038 
2001 


2000 


2000 
2001 


2100 
2100 


2100 


2000 
2001 
0254 
0255 
2100 


2100 


0422 
2004 


2100 


2000 


2004 


2001 


5000 


2001 


2002 


2004 


2000 
2001 
0564 
0567 
2000 
0371 


0570 


0256 
0257 
2002 
20085 
0400 


0440 


2004 
0410 
0426 
2004 
20035 
2001 
0415 


2000 


2004 
2001 
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.763 X 10° in 2000, 2001 


no ~0365, yos ~0367 


ba in 0236, 0237 


bz in 2002, 2003 


bz # bg in 2000, 2001 


pak 





0415 


0414 


0415 


0416 


0417 


0420 


0421 


0422 


0425 


0424 


0425 


0426 


0427 


04.30 


0431 


0452 


04.54 


04.55 


0436 


0457 


0440 


0441 


0442 


0445 


0444 


sl 


ad 


tm 


7] 


sf 


sf 


ta, 


ad 


ad 


ad 


ad 


35 


an 
30 


el 


o4 


02 


o4 


O02 


o4 


o4 
ol 
o4 
ol 
39 


o4 


oO 
o4 
55 


35 


2002 
3000 
2001 
2001 
2001 
0425 


3000 


3000 


2001 
2003 
5000 
2000 
9000 
2002 
2002 


5000 


2100 
2100 
9000 
2100 


2100 


2100 
2100 
0424 
0424 
0425 
2000 


2100 


2100 


2100 
2100 
2100 
2001 
2100 
20035 
2000 


2100 


0252 


2100 


02350 


0251 


2000 
0405 
2001 
2001 
2001 
2000 
0407 


0030 


0001 
Q001L 
0431 
0435 
0410 
2000 
0427 
2002 
0411 


04.05 


2002 
2005 
0400 
2002 


2005 


be in 2002, 2003 
bo #b3 #b4 in 2000, 2001 


} by in 2002, 2003 
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0445 
0446 
0447 


0450 


0451 
0452 
04953 
0494 
0455 
0456 
0457 
0460 


0461 


0462 
0463 
0464 
0465 
0500 
0501 
0502 
0503 
0504 
0505 
0506 
0507 


0510 


ad 


ad 


ad 


ta 


pr 
ht 


ad 


ad 


ad 


CX 


sf 


ad 


: 


09 


36 


34 


09 


ore 


od 


oe 


30 


300 


32 


o5- 


35 


39 


of 


OF 


O4 


30 


OZ 


Ol 


00 


20 


oF 


O2 


50 


3000 
2100 
2100 
3000 
2000 
0157 


2000 


2001 
0537 
2004 
0461 
2002 
0201 
0276 
3000 
3000 
0050 
2009 
2100 
2000 
2001 


3000 


O212 


2100 
0262 
0265 
2100 
0227 
2000 


0454 


2000 
O2ei 
0156 
O277 
2100 
0200 
0200 
2100 
2100 
0506 
0507 
[i] 
2002 
20085 


2100 


T7777 


0506 


0400 
2002 
2008 } 
0400 
2000 


0455 


0002 


2002 
2004 
0461 


0461 


0200 
0012 
1100 
0010 
20095 
0502 
2002 
2000 
2001 
0054 


0014 


2005 


48 


bl & bo- bz #b4 in 2000, 2001 
02545 in 2002, 2003 


t 
dws 


ap in 2000, 2001 


fo 


dvri 


(G41) * 2222 £41) x 27** ee 
iXQ-24 


test mx 1>i 


transferred to intecsration sub-routine 


U.T. to integration sub-routine 


[zi] in mg of 2005 


a, in 2000, 2001 





O511 
0512 
0515 
0514 
0515 
0516 
0517 
0520 
O21 
0522 
0525 
0524 
0525 
0526 
0527 
05350 


0531 


Ox 


sf 


ad 


Cox 


sf 


ad 


ex 


sf 


ad 


mr 


tm 


32 


Sl 
Sts) 


20 
o4 
50 
52 
ol 
59 
20 
54 
50 
o2 
ol 
55 
20 


o4 


2005 
2100 
2000 
2001 
3000 
0107 
2005 
2100 
2000 
2001 
3000 
0120 
2005 
2100 
2000 
2001 


5000 


0507 0512 
[ zi | 2002 
2002 2000 
2003 2001 
2100 0101 
0506 2005 
0507 0520 
[zi-2] 2002 
2002 2000 
2003 2001 
2100 O14 
0506 2005 
0507 0526 
[wa 2002 
2002 2000 
2003 2001 


2100 0125 


} ao in 2000, 2001 


\ az in 2000, 2001 


t az in 2000, 2001 


(cells 0532 through 0537 are filled with zeros) 


0540 
0541 
0542 
0543 
0544. 
0545 
0546 


0547 


sm. 


Cz 


sf 


ad 


50 


o2 


ol 


58 


29 


o4 


05 04 
2005 
2100 
2000 
2001 


5000 


0506 2005 
0507 0542 
ri A | 2002 
2002 2000 
2003 2002 


2100 0311 


by in 2000, 2001 


49 





0550 


0556 


0557 


Od561 
0562 


0565 


0565 
0566 
0567 
0570 
0571 
0572 
0575 
0574 
0575 
0576 


0577 


tm 


Ox 


sf 


ad 


Ox 


sf 


ad 


lH 


50 


Of 


ol 


35 


rae) 


o4 


30 


O2 


OL 


55 


25 


54 


0557 
2005 
2100 
2000 
2001 


5000 


0352 


2005 


2100 ws-1] 


2000 
2001 


3000 


0564 
2005 
2100 
2000 
2001 


5000 


0506 

0507 
’ 

[ wa | 

2002 

2005 


2100 


0506 


0507 


2002 
2005 


eL00 


2005 
0952 
2002 
2000 
2001 


0544 


2005 
0562 
2002 
2000 
2001 


0557 


2005 
0572 
2002 
2000 
2001 


0571 


\ bo in 2000, 2001 


} bz in 2000, 2001 


20 





1000 
1001 
1002 
1003 
1004 
1005 
1006 
1007 
1010 
1011 
1212 
1015 
1014 
1015 
1016 
HOLT 
1020 


1021 


Sl 


ad 


bl 


Ox 


CX 


ox 


sl 


pr 


ra | 
39 
35 
05 
Oe 
20 
O2 
O2 
ral 
o4 


21 


o4 


O02 


02 


Ol 


0000 
1014 
2000 
0000 
2002 
2002 
2002 
2100 
200] 
2001 


2000 


5000 


5000 


f 


7400 


1016 
2001 
2100 
3000 
L000 
2007 
2001 
2001 
2005 
2100 


1020 


2100 


f 


f 


2001 
10135 
1022 
1014 
2000 
2002 
2000 
2002 
2001 
1005 


0001 


0004 

00350 
a 
i 


f 


(cells 1022 through 1047 are filled with zeros) 


1050 


1051 


ad 


ad 


1052 [ ad 


1053 


1054 


1055 


1056 


tm 


ad 


oS) 


ere) 


O90 


o4 


ore) 


o4 


O4 


2100 2100 
1060 1057 
1305 Ai 2100 
3000 2100 
1062 1057 
1061 1057 
3000 2100 


1057 clear i tally 

1052 

2000] 

1000 Uet. to convert and print sub-routine 
1057 ifloi 

1051 


1160 


ol 





1057 


1060 ad 


1061 


1062 


bo 


1505 


0050 


0001 


2100 


if 


€ 


2000 


nf 


ue 


i X grie 


mel ff 


(The following is the sub-routine for the integration of n lst order 


equations by the Runge-Kutta-Gill method) 


1100 
1101 
1102 
1103 
1104 
1105 
1106 
1107 
1110 
1111 


11i2 


sl 
ad 
tm 
tm 
ad 
ad 
ad 
ad 
ad 


[sm 


ad 


1115 [ su 


1114 


ad 


1115 | mr 


1116 
ee 
1120 
1121 
1122 


1123 


8u 


to 


ad 


ad 


ot 
55 
o4 
o4 


59 


5D 
55 
59 
50 
59 
56 
59 
25 
o6 
o7 
20 
oO 
59 


30 


5000 
1246 
11635 
5000 
2100 
2100 
2100 
1250 


1254 


2004 
2004. 
2006 


2006 


1251 


1251 


1272 


1242 
2000 
1303 
2100 
2100 
2100 
1177 
1251 
1252 
Bd 

125) 
2001 
1250 
Ad 

2005 
2100 
1301 
1242 
1270 


NZL 


2000 
1160 
1200 
1107 
1247 
1250 
1251 
1252 


gp eva 


2001 | 


LS 


2007 
1270 
Agel 


1124 


Prepare exit to main program 


Test if first time through 
Teeter 
O-~+k (clear print tally) 
O +j (clear iteration tally) 
1 xX 27) (rosct var. noe tally) 
i xX 2-12 4 5 x 2784 tally 
G Jj in M, ana [oJ in M, of 1111 
Bx q: 
[2] in My of 1113 
1-84) 
[A'] in ig of 1135 
ACG-8 qi) 
(In )(x}- 893 = = AX} /z 
test | r|>] 
(ANC Big?) Te AX; 
i xX 2795 from i x 27 l? 
ix arte of i x 2796 


[ x0] in li, and M, of 1124 


(truncation test) 


02 





1124 
1125 
1126 


i 


ad 


Sm 


ad 


ad 


1130 [su 


od 
1132 
1133 
1154 
1155 
1156 
eo 
1140 
1141 
1142 
1143 
1144 
1145 
1146 
1147 
1150 
fo) 
LioZ 
1153 
1154 
1155 


1156 


ad 
E 
sm. 
sm 
ad 
mr 
ad 
sm. 
ad 
[sm 
ad 
tn 
ad 
ts 
ad 
tn 
ad 
ts 
ad 
tm 
ad 


ts 


59 
30 
09 
re) 
o6 
39 
20 
30 
30 
39 
20 
39 
30 
ee) 
30 
39 
04 
39 
iy 
ore) 
O4 
O90 
ive 


O5 


Ste) 


17 


1251 
L251 


1257 


1241 
sed) 

2007 
2000 
2000 
2006 
2001 
1111 
1276 
2002 
1177 
1300 
2100 
2010 
1250 
1176 
2100 
2020 
1163 
1302 
2100 


2040 


2007 
1255 


20035 


L274 
1274 
2007 
1275 
2006 
1242 
2003 
BATE 
1251 
l2o1 
ae a 
3000 
1175 
1250 
2100 
3000 
1247 
1247 
2100 


3000 


J 
xi J 


2005 


1273 


11350 


2007| 


1132 


2000 | 


2007 
2000 
2001 
2006 
2002 
2005 
1142 
7. | 
P29. 
1107 
1251 
1164 
1250 
1160 
1250 
1166 
1247 
L160 
1247 


1170 


xd TA] xj Al ; xd 
i a ie 1 
i X 2724 from i X Q-le2 
iX gle / i x 2724 
iy, JT 
£43] in My and [YJ in My of 1130 
é Ae 
4 ’ 
q,~ A; 


2 
[ x; J in My and Ga in My of 1132 


Xo x C 


iy) x enh 
(4, z;) 
ieme. X 8 

id i P| 
Laer (c’-!) %; ] X8 


37; x ent 


fale Coit) af ]xB 4 307 ABs 1 xg" 

ee) 

if in Mz of 1142 
Laid { al 3 
(GM Ka a 4; 


in Mz of 2003 


(441) X orl2 + 4 x grle 
test (n #1) >i 

1x 19712 + 4 x 272 

print j X 27°* if ts 2010 up 
Gime 4 X 2-6 
test 4 > j 

O > 4 x 2744 

accumulate © if 2020 up 
k 4 ees ie 

test K>k 

O—-k 


print X; octal if 2040 up 


55 





1157 
1160 
1161 
1162 
1165 
1164 
1165 
1166 
1167 
1170 
nel 
ue yg 
1173 
1174 
WOT 
1176 
ET 7 
1200 
1201 
eee 
12035 
1204 
1205 
1206 
1207 


1210 


tra 


pr 
ad 


ts 


[= 


ad 
ad 
ad 


ad 


ad 
ad 


ad 


04 


Od 
54 
Sts) 
54 
uy, 
ol 
54 
oi 


o4 


55 
5O 


59 


50 
59 
05 
38 


30 


5000 


1250 
3000 
1501 
3000 
2020 
1505 
5000 
1504 


3000 


0001 
1250 
i208 
1500 
1255 
2000 
2002 
1256 
1240 


2002 


2100 


2100 
2100 
1504 
2100 
3000 
2100 
2100 
2100 
2100 
0001 


0004 


Hy 


1500 
2000 
2000 
2002 
1235 
2005 
2004 
2000 


1242 


1050 


0001 
0001 
1147 
1504 
1153 


1175 


1160 


n fi | 


1160 


2000 
ee, 
2002 
1234 
2003 
2004 


1257 


U. T. to new print routine 


U. Te to X, SR 


1 x 9736 

print j X 274 

ieee 

Xo (mA1)T — XA mbt 
Ucmeie 


Bit ge Seo cs ssi iy, 


print A, = t, Xy, Xo, oeeeX, 


1 X 2744 
ye phe 

1x 2712 

[X,] in ty of 2000 

[xX,] im My of mod. dumy 1 

(9, ] im ty of 2002 

[4,] in M, of mod. dummy 2 

[z,] im Mg of 2008 

[ q.] in My and [X.] in Me of 2004 
Log) Se Bh ees [% Jin U, of mod.dumy 3 
CX,] im My of mod. dumy 4 

[9,] in Mg of 2002 


o4 





1211 sm 30 1300 1242 2000 eee iy 2000 

1212 ad 35 1245 2000 1173 sets up print of X,, 1, eee. X, 
1213 su 36 2000 1163 2001 n-X 27°8 in 2002 

1214 ad 35 1244 1+\2001 1171 sets up print of X,,X5 , ees. X, 
1215 ad 35 2100 2100 1270 0 i x 9756 

1216 ad 35 2002 1270 2008 [4:] in lg of 2005 

1217 ad 35 1245 2005 1220 [4,] in Mg of 1220 

1220 [ad 35 2100 2100 q; J] oa 

1221 ad 35 1163 1270 1270 (441) X 2736 —» 4 x 2756 
1222 tm 34 1270 2000 1224 test i pn Al 

1223 tm 34 3000 2100 1216 

1224 ad 35 2100 1163 1303 signal for first time through 
1225 tm 34 3000 2100 1104 apy 

1226 pr 21 2006 1227 0001 print _ 

1227 ht 22 fc f f halt 

1230 1304 isk f [x] x 272? 

1251 su 36 £ 2001 2004 dumny 1 

1232(su 36 X, = 2002S 2004] = modified dummy 1 for 1112 
1233 sm 30 f 1253 2001 dummy 2 | 
1234 [sm 30 gq, 1253 2001 | modified dumy 2 for 1111 
1235 02 f f 0014 

1256 su 36 if i 2007 dummy 3 

1237 [su 36 qo X, 2007 / modified dumy 3 for 1130 
1240 mr 26 f 1263 2000 dummy 4 

1241 [ m 25 X 1263 2000] modified dumy 4 

1242 02 f f 0030 
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1243 
1244 
1245 
1246 
1247 
1250 
1251 
1252 
12535 
1254 
1255 
1256 
1257 
1260 
L261 
icon 
1265 
1264 
1265 
1266 
1267 
1270 
1271 
1272 
1275 
L274 


L275 


pr 


pr 
ad 


ad 


al 
21 
ore) 


o4 


O02 


02 


25 


ors) 


02 


1504 
1505 
2100 


3000 


4000 
5520 
9920 
6525 
4000 
5520 
5520 
4000 


2004 


1504 


3000 


2100 
2100 
2100 


2100 


2007 
2007 


2020 


2007 


2007 


1257 


2007 


OOOL 


OOOL 


4151 


4151 


2929 


4151 


4151 


2005 


1504 


0005 


f 


dummy print for 1175 
dummy print for 1171 
dumny for 1220 
dummy exit coamnand 

k X 2786 (print tally) 

j X 27° (iteration tally) 
i X 27-12 (var. no. tally) 
ix2l2 y 5 x 2% 

0 

Be 

Be 

Be 

A° = 1/2 

als AZ 

aes AP 

A° = 5/6 

cos 1/2 

os 2 

cfs - VAP 

co = 1/2 

dumy for 1115 

i ee? 

Lat gee OSE rae 
dumny for 1124 


ix 2t24 i x 2°24 


3 x gvt 





1276 sm 50 2002 1277 £ 

1277 f f 0003 

(Cells 1300 through 1332 were filled by tape. These cells contain the 
values of t, 23, Wy, etoe These values are not included here because they 
depend upon the values of computed parameters--see Chapter V.) 

Some notes on the program: 

All addresses appearing in brackets are not put in by punch cards but 
rather by other parts of the program. They are included for purposes of 
clarity in deciphering the programe 
The following notes apply to the intesration suberoutine: 

ts 2010 -- up for iteration number tally 

ts 2020 -- up for goneration of x, 
ts 2040 -- up for print of x; octal 
n * no. of dependent variables (equations ) 


(i 


K 


incremental change of independent variable 


readout at every K point 

X,# independent variable 

New print or readout routine to start at 1050. 

Automatic print of cell 2006 indicates interval of integration too 
large. 


Error of integration is proportional to the 5th power of T. 


o7 








CHAPTAR V 


CONFLUDING REMARKS AND RECOLIMENDAT IONS 
FOR PURTHER WORK 


As noted previously, it would be fatuous at this point to draw any 
definite conclusions anent tho solution to our problem proposed in the 
first chapters On the basis of the computer solutions we cannot say with 
absolute certainty that the results are compatible with experimental 
observations. We can say, howevor, on the basis of our sketchy results 
that the theory and method are promising. Without equivocation we can 
say that the computer program is eminontly successful--it does work and it 
does solve the equations proposed as a description of the traveling 
density wavese A great deal of work is forecast in ordor to establish 
the identity of these equations with the phenomenon in the positive column. 
Great numbers of initial conditions should be inserted in the program and 
the coefficients should be varied in order to eventually establish the 
validity of the equations. It is bolieved that if this is done, the 
resulting solutions by the computer will correspond to the observations 
made under actual conditions in the laboratory. lIlowever, no precise 
conclusions can be drawn fram a very few computations. These computations 
cannot be made in a few hours. It is prodicted that the collection of 
enough data will require on the order of several hundred hours of computer 
time. The important point to be remombered is that a workable progran is 
availablo and patience and persevorance will reward the eventual researcher 
with the answere It appears from initial computer runs that the term 
zng is dominated in the equations (4) by other terms when pts Noe 


The initial computer runs suggest that "floating point arithmetic" might 
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be necessary throughout the program due to the drastic changes noted in 
derivatives and the time required for computation using the present 


programe lio other trends were noted from the initial oomputor runs. 
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Varney, Re Ne 


Wornbeck, Je Ae 


Chanin, Le. Iie and 


Biondi, lie Ae 
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Farris, Ve De 


Grahan, We Je and 
Ruble, Ae aie 
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